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We investigate the transport properties of a classical wave propagating through 
a quasi-periodic Fibonacci array of waveguide segments in the form of loops. The 
formulation is general, and applicable for electromagnetic or acoustic waves through 
such structures. We examine the conditions for resonant transmission in a Fibonacci 
waveguide structure. The local positional correlation between the loops are found 
to be responsible for the resonance. We also show that, depending on the number 
of segments attached to a particular loop, the intensity at the nodes displays a 
perfectly periodic or a self-similar pattern. The former pattern corresponds to a 
perfectly extended mode of propagation, which is to be contrasted to the electron or 
phonon characteristics of a pure one dimensional Fibonacci quasi-crystal. 

I. INTRODUCTION 

Photonic band gap (PBG) materials have remained an active area of research in condensed 
matter physics and materials science over the past several years With 
these systems one comes across the interesting possibility of influencing the propagation of 
electromagnetic wave by creating gaps in the band structure of synthetic periodic dielectric 
structures jj] . Such studies constitute an important part of mesoscopic physics, which focus 
on the possibility of localization of light P, Si , a direct evidence of which has already been 



reported for a strongly scattering media of semiconductor powders 



Networks formed by segments of one dimensional waveguides provide examples of an 
alternative class of PBG systems which do not require a material with large dielectric con- 



stant. Zhang et al 



have experimentally observed Anderson localization of light in a three 



dimensional network, composed of coaxial cables which played the role of a practically one 
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dimensional waveguide. Vasseur et al 19] have investigated the photonic band structure of 
a comb-hke waveguide geometry composed of danghng side branches grafted periodically 
along a mono-mode waveguide. Their work represents a one dimensional photonic crystal 
enabling one to investigate the occurrence of localized states in such systems. The network 
system is capable of producing large gaps even in one dimension, though it may be sensitive 
to the structure of the unit cell. 



context of propagation of electromagnetic and acoustic waves |lOllllj. Experiment has been 
done [l^ on systems consisting of loop structure made by slender tubes which are pasted 
on similar tubes of finite length. These circuits allow the propagation of electromagnetic or 
acoustic wave through practically a single channel. A theory for one dimensional network 
thus becomes relevant. However, the results presented so far mainly deal with periodic 
systems, and the influence of a variation from periodicity in the geometry of the network on 
the transport of classical waves practically remains unaddressed. 

Motivated by this observation, we undertake a systematic investigation of the transmis- 
sion of classical waves through a q uasi- p er io dicallv ordered Fibonacci SLS. We use a model 
proposed initially by Zhang and Sheng [1^ to examine wave propagation in a random array 
of waveguide segments in the form of loops. Two consecutive loops join each other at a 
node. Scattering occurs only at the nodes, each pair of which holds an arbitrary number 
of segments between them. A one dimensional Fibonacci lattice with the golden mean irra- 
tionality is already known to exhibit a purely singular continuous spectrum, free from the 
existence of any extended wave function The interplay of quasi-periodic order and the 
number of waveguide segments between the consecutive nodes (each segment providing a 
propagating channel) in a closed loop geometry is likely to produce interesting transmis- 
sion characteristics, which is the prime objective of the present study. Secondly, whether 

extended modes of classical wave propagation in a quasi-periodic geometry exists, is yet to 

j 

be known, to the best of our knowledge. Such extended states do exist ^] for an electron 
traveling through a certain group of quasi-periodic lattices as well as in random lattices 



Recently, network models based on serial loop structures 
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with correlated disorder where one identifies local clusters of atomic sites, causing 

resonance and extended wave-function. However, delocalized eigenmodes in the context of 
quantum transport does not occur in a golden mean Fibonacci chain. It is therefore, in- 
teresting to see whether, multiple loops have any role to play in this regard, in the case of 
classical wave propagation. 

We find that the band structure of an array of Fibonacci SLS has a strong dependence 
on the number of waveguide segments between consecutive nodes, which is refiected in the 
transmission coefficient across a finite Fibonacci SLS. We have also been able to identify 
local clusters in a Fibonacci SLS, which are responsible for producing resonance leading 
to finite transmission. The transmission coefficient, at special values of the wave- vector 
exhibits cyclic variation as a function of the system size. The intensity distribution at the 
nodes corresponding to such cases may exhibit a completely periodic pattern, or a self- 
similar (fractal) pattern depending on the mutual relationship between the lengths and 
the numbers of the waveguide segments in the two basic loop structures constituting the 
Fibonacci geometry. 

In what follows, we describe our model, method and the results. A comparison with the 
corresponding periodic system is made wherever necessary. The formulation remains valid 
whether we have an electromagnetic wave or sound wave propagating through the network. 



II. THE MODEL AND THE METHOD 



Following Zhang and Sheng we consider a network formed by waveguide segments of 
fiexible length joining together to form a Fibonacci sequence of loops A and B [Fig. 1(a)]. A 
binary Fibonacci sequence is generated from a seed A, by using the infiation rule, A AB 

n 

and B ^ A [13]. The 'wave-function' i^i^i+i within any segment of length U^i+i between the 



nodes i and i + 1 is given by 



sin[fc(/i,i+i - x)] sinA;a; 

= sinH ■ . + ^-+1 inH.. , 

sm /ctj^j_|_i sm /ctj^j+i 
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where, ipi is the amphtude at the ith mode, and k is the absolute value of the wave vector 
in the loop. For example, for an electromagnetic wave, k = ^^^^ where ui and Cq are the 
frequency and the speed (in vacuum) of the electromagnetic wave respectively, e is the 
relative permittivity of the dielectric medium, which we may assume to be real. A complex 
dielectric constant can easily be dealt with. The continuity of the wave function at the nodes 
and the flux conservation criterion are used to map the problem of wave propagation in such 
loops into an equivalent tight binding problem of electron propagation on a one dimension 



lattice 



12| . The resulting difference equation is, 

(E - ei)^i = ti^i+itpi+i + tj_i,i'0i_i (2) 



where, we parametrize E = 2 cos kli^i^i and, define the effective on-site potential 

Q = 2cos(A;/,,+0 + J2 cot^S, + Yl cot^Si (3) 

m=l m=l 

The effective nearest neighbour 'hopping integral' is given by, 

Ni.i+i ^ 

^^'^+1 = E . M (4) 

m=l 

Here, Ni^i±i is the number of segments between the nodes i and i ± 1. O^^^l-^ = kl^^^^]^^ is 
the corresponding phase acquired in the w}^ segment. The problem of wave propagation 
now becomes mathematically equivalent to the problem of transmission of an electron with 
energy E through the mapped lattice, where, —2<E<2. The set of equations (2) is 
conveniently cast into a matrix form, 

MA \ (5) 




where, 



is called the 'transfer matrix' 



1 



In a Fibonacci SLS of two loops A and S, the loops may differ from each other either 
in terms of the lengths of the segments, or, in the number of segments, or in both. Using 
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the language of electron propagation, on the mapped lattice, which now becomes a one 
dimensional Fibonacci chain, the 'on-site potentials' jl^ assume three distinct values (in 
respect of the nearest neighbour loops), viz, 

= 2 cos 9i + m[cot Oi + cot 62] 

e/3 = 2 cos 6*1 + m[cot 9i + cot 62] + n[cot 9^ + cot ^4] 

^7 = e/3 (6) 

There are two kinds of 'hopping integrals' II and ts for the equivalent electron problem 

which are arranged accordingly to the Fibonacci sequence. These are given by, 

m m 

tL = —^ + 



sin 6*1 sin 62 

IT' IT' 

ts = —pr + 7 
sm 6^ sm 6/4 

Here, we have assumed that the A-loop has m segments of length ai in the 'upper' branch, 
and m segments of length 02 in the 'lower' branch. The S-loop has n segments in both 
branches, the lengths being as and 04 for the upper and the lower branches respectively, and 
6j = kaj with Cij — CLi, 0,2, CL3, CLa- The asymmetrical Fibonacci SLS can thus be easily studied. 
On the equivalent one-dimensional lattice [Fig. 1(b)], we define three transfer matrices, viz 











M„ j 






ts ts 




I 1 i 


I 1 



E-ej _ts_ 
tL tL 



1 

corresponding to the sites a, (3 and 7 respectively fl^. We have examined the spectral 
properties of both the infinite and the finite Fibonacci SLS for various combinations of m, 
n and a/s. The fundamental features that we propose to investigate are revealed even if we 
keep all the lengths identical. Therefore, without losing any generality we stick to the case 
where ai = a2 = a^, = = a and present the results in the following section. 

III. RESULTS AND DISCUSSION 
A. Band Structure 



The first few generation in a Fibonacci sequence of A and B are as follows: 5*1 = A; 
S2 = AB; S3 = ABA; = ABAAB and so on. An infinite Fibonacci sequence of loops 
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FIG. 1: (a) A portion of an infinite Fibonacci array of waveguide loops, (b) The equivalent Id 
Fibonacci chain. The smallest resonating clusters (3"f and a have been marked. 

A and B (or oi and ts) may be recursively generated by periodically repeating rational 
approximants of increasing size. By noting the individual sites (in terms of the equivalent 
one-dimensional chain), one can calculate the product transfer matrix across an genera- 
tion Fibonacci rational approximant. The condition for a certain k-value to be an "allowed" 
one is, given by | |< 1, where, xi = \TrMi Mi being the transfer matrix for the 
jth generation Fibonacci lattice, and satisfies the recursion relation Mi = Mi_2Mi_i with 
Ml = Ma and M2 = M^jj = M^Mjj. The relation between x^'s is, 

Xi = Xi-2Xl^l - Xl-3, I > 4 

This is called the 'trace map'. Using the trace- map method we have studied the distri- 
bution of the allowed values of the wave- vector A; as a function of Am = m — n keeping the 
lengths of all the segments in each of the loops A and B equal. For this we have selected 
an arbitrary rational approximant of a Fibonacci lattice, which consists of 55 loops. The 
case Am = corresponds to the perfectly periodic arrangement of identical loops, for which 
we have a continuous distribution of "allowed" eigenvalues from A; = to 27r (as we have 
shown) with a = 1 in arbitrary units. With n = 1 [Fig. 2 (a)] the spectrum starts splitting 
into three major sub-bands, with a couple of minor clusters of allowed k-values peeping 
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in between the major sub-bands. These clusters do not grow when Am exceeds a certain 
value. Each major sub-band gets split up, exhibiting a fragmented structure as m increases 
keeping n fixed at 1. However, when Am becomes quite large, the further fragmentation of 
the major sub-bands becomes hard to detect numerically. The features are sensitive to the 
starting value of n. Fig.2(b) shows the band structure, where the number of segments in 
the B-\oop is 2n = 10. Here, the minor (weak) sub-clusters around k = 1.5 and k close to 
5 grow in width, as well as in length. With increasing values of n the spectrum assumes a 
very prominent four-sub-band structure. The self-similar splitting of the spectrum is there 
as Am increases, which of course is limited by the resolution of the machine as Am increases 
to large enough values. 



B. Transmission coefficient 

For computation of the transmission coefficient, we adopt the widely used technique ^| 
of fixing semi-infinite leads to the extreme nodes of the equivalent one-dimensional Fibonacci 
lattice obtained by discretizing the wave equation. The lead is described by identical and 
equispaced sites with on-site potential eo, set equal to zero and connected to each other via 
nearest neighbour hopping integral to, which we set equal to unity. An equivalent problem of 
an electron with energy E = 2 cos ka (in unit of to) travelling along such a lead and entering 
the finite sized sample will give us the transmittance of the actual problem of classical waves 
propagating through the Fibonacci SLS. Following Ref. the transmission coefficient of 
a finite equivalent one-dimensional Fibonacci segment consisting of sites a, P and 7 and 
hopping integrals ti, ts is given by, 

~ [(A2 - P21 + (Pu - P22) cos ey + (Pn + ^22)2 sin^ 9] ^ ^ 

where P = Yli M^. Mj are the transfer matrices corresponding to each site (including the 
boundary sites) in the equivalent chain that is clamped between the leads. Here 6 = ka. We 
have used the formulation first to get an idea of wave propagation through a periodic array 
of identical loops, each consisting of the same 'm' number of segments of equal length. 
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The transmission coefficient for such loops in periodic arrangement turns out to be, 

T — (9) 

[2m(C/Ar-2(x) - ?7^(x)) + ^i=:|^?7^_i(x) sin2^]2 + (1 + 4m2)2i72 _,(^^ 

where, Un{x) is the A^th order Chebyshev polynomial of second kind, a = ai = a2 = = 
and X — — cos ka — — cos 6. 

In Fig. 3 we plot the transmission coefficient for periodic array of 55 identical loops. It 
is interesting to see that in the present model the loops touch each other, and there is no 

real gap opening up in the spectrum, except for the isolated zeroes at ka — 0, tt, 27r, The 

envelope of the values of T between consecutive zeroes can be worked out from (9) to drop 
following a power law, viz, T ~ ^ with increasing values of m, when m is large. 

The calculation of transmission coefficient for a Fibonacci array of two different loops 
brings out an interesting difference. Fig. 4 shows the transmission coefficient against the 
wave-vector for a 9th generation Fibonacci sequence of 55 loops. The loops differ only in 
the values of the number of waveguide segments m and n, n being fixed at 1. We have kept 
Oi = 02 = as = 04 = a as before, and results are shown for < fca < 27r with m — 2 
[Fig.4(a)]. The corresponding lengths of the wave-guide segments with frequency range for 
an electromagnetic wave for example, can easily be extracted to design a possible network 
for the light localization in a quasi-periodic network. The gaps in the spectrum are quite 
clear. With m taking up values 4 [Fig.4(b)] and 6 [Fig. 4(c)], the gaps clearly widen and the 
transmission spectrum exhibits a global three sub-band structure. It should be appreciated 
that the transmission spectrum is sensitive to the absolute values of the parameters of the 
system. For example, in any given generation, by increasing the lengths of the segments for 
any loop A, or or by increasing the value of n for a fixed m (> n) one finds the growth 
of new sub-bands in the transmission characteristics. However, the additional sub-bands 
disappear when one increases the difference between m and n appreciably. 



For certain values of the wave- vector, appropriate choice of the lengths of the segments 
may lead to interesting fluctuations in the transmission spectrum. This behaviour is, how- 
ever, sensitive to the number of A and B loops in the Fibonacci array. To clarify, let us 
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define 9j = kaj, j =1,2,3 and 4. It is simple to work out that, for 6'i = | and ^3 = ^, witfi 
Oi + 02 = {2p + 1)11 and 63 + 64^ = {2q + l)7r, p and q being integers, eacli of the transfer 
matrices M^, Mg and becomes equal to — zcr^, provided one adjusts m and n such that, 

msin^s = nsin^i. (10) 

In the equivalent electron problem we now have E = 2 cos^i = and = with ti = ts- It 
implies that the equivalent one- dimensional chain becomes indistinguishable from a perfectly 
ordered chain of pseudo-atoms. In this situation, a straightforward calculation shows that, 
a Fibonacci array of loops at the 31 and (3/ + 1)*'^ generation, with / > 1, will exhibit a ^ 
decay in the transmission coefficient. For all other generations, with the above specification 
of parameters, we get T = 1. 



C. Dimer like correlations and resonance 



Looking at the equivalent one dimensional Fibonacci chain (Fig. lb), we find that the 
cluster f3'~f (put in a box) appears in pairs, as well as in isolation. The site a (put in a circle) 
however, is always single. The /?7 — f3'~f pairs, whenever they appear, are always fianked by 
two a sites. Now, if the traces of the 2x2 transfer matrices M^^ and Ma can be made 
to vanish simultaneously for a certain combination of the system-parameters and the wave 
vector, then the cluster a — — 7 — /9 — 7 — a will contribute an identity matrix to the 
entire product of the transfer matrices at any given generation of the Fibonacci chain jisl ]. 
This leads to a resonance in the above cluster. As has been explained elsewhere 3|; such 
resonance then takes place locally throughout the chain, reducing the long product of transfer 
matrices, viz, M^+g to Mi for / > 1. This also explains the six cyclic invariance of the transfer 
matrices in a Fibonacci chain However, the sizes of the minimal clusters responsible for 
resonance depend on the choice of the wave vector, and definitely, on the combination of the 
other parameters of the system. To clarify, we work out a specific example. Let us select, 
61 + 02 = (2p + l)7r and 63 + 9^ = {21 + l)7r, p and q being integers, equal or unequal. This 
choice automatically ensures ea = ep = e^/ = E and a six-cycle of the transfer matrices for all 
/c- values consistent with the above conditions for any suitable combination of the a^'s jlsl ]. 
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In particular, when the special condition (10) is satisfied, it is simple to check that ti = ts 
and Ma = Mp = = —iay. The products M^^ = = —I, and the resonance condition, 
as discussed above, is trivially satisfied. With (3'y — P'j and a — a playing the roles of the 
minimal clusters responsible for resonance. Otherwise, if m sin ^ nsin^i i.e. 7^ ts, we 
get back the well known pure transfer model of a Fibonacci chain jl^j with E = 2 cos 9i at its 
band-centre. Again, a six-cycle wave-function results, and the intensity | "V^i P has a purely 
multi-fractal character • The multifractality can be removed by controlling the values of 
m and n. For example, with kai = | and ka^ = ^, we get ti = ts^- So, m = 2n makes the 
equivalent one- dimensional lattice completely periodic. \ ipi P is also periodic in this case, 
though the loops are arranged in a Fibonacci sequence. We get an extended eigenmode. For 
all other cases with m = 2n, the transfer model, and hence the fractal nature of | ■j/'j p at 
the nodes is restored. This feature is illustrated in Fig. 5. 



IV. CONCLUSION 



We investigate the propagation of classical waves in a Fibonacci array of waveguide 
network. For specific considerations of the system parameters one can have extended, non- 
fractal wave function is the system, contrary to the quantum counterpart of the problem. 
Resonant transmission in such network can be attributed dimer-like correlation in the place- 
ment of the loops. The study shows that one can design such a network as an example of 
quasi-periodic photonic band gap system. 
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